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Abstract. Recent analyses of the cosmic microwave background (CMB) maps 
from the WMAP satellite have uncovered evidence for a hemispherical power 
anomaly, i.e. a dipole modulation of the CMB power spectrum at large angular 
scales with an amplitude of ±14 percent. Erickcek et al have put forward an 
inflationary model to explain this anomaly. Their scenario is a variation on the 
curvaton scenario in which the curvaton possesses a large-scale spatial gradient 
that modulates the amplitude of CMB fluctuations. We show that this scenario 
would also lead to a spatial gradient in the amplitude of perturbations and 
hence to a dipole asymmetry in any highly biased tracer of the underlying density 
field. Using the high-redshift quasars from the Sloan Digital Sky Survey, we 
find an upper limit on such a gradient of |V<rg|/cr8 < Q.027t\~J (99% posterior 
probability), where ri ss is the comoving distance to the last-scattering surface. 
This rules out the simplest version of the curvaton spatial gradient scenario. 



1. Introduction 

The simplest inflationary models for the origin of structure of the Universe has made 
a number of successful predictions, including the flatness of the Universe {VLk = 0), 
and near-scale invariance (n s w 1), adiabaticity, and Gaussianity of the primordial 
fluctuations. Another prediction of these models is that the perturbations should be 
statistically homogeneous and isotropic. It was therefore a surprise when, after the 
cosmic microwave background (CMB) maps from the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite were released [U HI several groups reported evidence 
for large-scale anomalies in the data that were suggestive of inconsistencies with the 
paradigm of a statistically isotropic Gaussian random field, or with the theoretical 
power spectrum. One widely discussed anomaly was an alignment of some of the 
low multipoles with a so-called Axis of Evil [4j [5l [6l [7]. Another was the CMB 
quadrupole C2 , which was unexpectedly low for the standard cosmological model [I] . 
A third anomaly was a difference in the observed CMB power spectrum between two 
hemispheres [8j OH [TQl [Til HI] • A variety of theoretical explanations for these anomalies 
have been proposed [El[H[T51[TJ[TIl[T8l[T9l[20]. 

The hemispherical power anomaly is the main subject of this paper. Its properties, 
as measured by Hoftuft et al |12j . can be summarized as follows: (1) The direction 
of greatest power is p = (l,b) = (224°, -22°) ± 24°, near the South Ecliptic Pole. 
(2) The amplitude of large-scale fluctuations (I < 64) is A = 0.072 ± 0.022 (3.3fi), 
meaning that the amplitude of the CMB perturbations At is 1.07 times the mean 



Hemispherical anomalies 



2 



value in the p direction, and 0.93 times the mean value in the — p direction^ (3) 
The asymmetry is present in at least three of the WMAP bandpasses (41, 61, and 94 
GHz), and in the foreground-suppressed internal linear combination (ILC) map [21J. 
(4) Only limited information is available about the ^-dependence of the asymmetry. 
There are indications that it continues at least to I — 80 [12] and perhaps beyond [11] , 
but other analyses have found no evidence for an amplitude asymmetry at the first 
acoustic peak [28] . 

The large-scale anomalies are a difficult subject, in part because there is no 
clean way to assess their true significance. While the quoted significances sometimes 
exceed 99%, some of the analyses are a posteriori. In order to make progress it is 
necessary to consider the theoretical explanations for the large-scale anomalies, and 
ask what new predictions these models yield [22]. For some of the proposed ideas 
that only affect the low multipoles, the CMB temperature field has reached its cosmic 
variance limit and further tests must await future CMB polarization data with lower 
noise and foreground residuals than WMAP. On the other hand, one class of models 
intended to explain the hemispheric power asymmetry make a host of new predictions 
testable with near-future or even current data. Erickcek et al [20] recently proposed 
a variation on the curvaton model (23] [211 Ell EE] : an inflaton 4>, which dominates 
the energy density during inflation, and a curvaton <j, which is effectively massless 
during inflation but decays after inflation ends to produce an energy density oc cr 2 [§| 
In the variation of Erickcek et al [20] . there was initially a large-scale gradient in a; 
since the primordial density perturbations are proportional to S(o~ 2 ) ~ 2aSa (assuming 
they are sourced by the curvaton; inflaton-induced density perturbations may also be 
significant), the amplitude of the density perturbations is statistically inhomogeneous, 
being proportional to the local value of A uniform gradient in the background field 
a then results in a dipole modulation of the CMB power spectrum in the direction of 
some unit vector p || V|<r|, with an amplitude given by the fractional variation of a 
across our horizon volume. 

As noted by Erickcek et al ED], the proposed model could be tested with future 
CMB data, especially with the upcoming Planck satellite, which will observe many 
more modes in the CMB than WMAP and should definitively confirm or reject the 
power spectrum asymmetry possibly seen by WMAP. (The model could possibly even 
be constrained further by using the higher multipoles in WMAP.) However such a 
model also makes predictions for large-scale structure: since the dipole asymmetry 
in the power spectrum varies slowly with wavenumber k, it generates a large-scale 
gradient in the amplitude of fluctuations ag. Thus the abundance of rare massive 
haloes should contain a spatial gradient. Thus the angular distribution of these 
massive haloes (or any luminous tracer thereof) should contain a dipole in the same 
direction p singled out by the CMB. Moreover, using the measured dipole anisotropy 
of the CMB power spectrum and our knowledge of the halo mass function, one can 
specifically predict the amplitude of these features as a function of redshift and halo 
mass. Such features, if detected with the correct dependences and alignment of the 
CMB, would be a smoking gun for new physics, whereas nondetection of the large 
scale structure dipole would rule out the simplest version of the Erickcek et al [20] 

X The asymmetry in the C;s is twice this, ±14 percent, since the power spectrum is the square of 
the amplitude. 

§ We have assumed here that the curvaton never contributes more than a small fraction of the energy 
density of the Universe; this is true in the Erickcek et al model, but it is not true in all curvaton 
models. 
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model. 

The purpose of this paper is two- fold: first, to develop the formalism to test for 
a spatial gradient in cr 8 ; and second, to apply this formalism to the current large- 
scale structure data and measure Vers- Our primary motivation was the Erickcek et 
al [20] model, but we note that any early-Universe explanation for the hemispherical 
asymmetry is likely to produce Vers and hence our measurement should provide tight 
constraints. [We note that a statistically anisotropic power spectrum P(k) cannot 
produce the dipole signature seen by Refs. [3 [TUJ [HJ [12] because it is invariant under 
k — ► — k.] For our measurement of Vas, we use the high-redshift spectroscopic quasars 
from the Sloan Digital Sky Survey (SDSS) as described below. 



2. Theory 

2.1. The model 

The Erickcek et al [20] model produces a local primordial fluctuation amplitude 
A-ji(k) that has a spatial gradient: 

A n (k,r) = A K (k,0)(l+p- — ), (1) 

V r lss / 

where r\ ss is the distance to the surface of last scattering. This normalization was 
chosen so that in observations of the CMB the amplitude of variations is modulated 
by a factor 1 + A cos (9, where 6 is the angle between the gradient direction p and the 
line of sight h and A = \p\ (with minor modifications due to reionization; see below). 
Constraints in this paper will be reported on the amplitude A and direction p of this 
modulation. 



2.2. Effect on CMB 

A large-scale gradient in the primordial fluctuation amplitude will produce a 
corresponding dipole in the CMB anisotropies. On small scales, where the added 
power from the integrated Sachs- Wolfe (ISW) effect and reionization are negligible, 
the local CMB power spectrum in a particular region of the sky n is 

C\ oc {h) = Cie- 2AT{n) (l+p-h) 2 , (2) 

where Ar(n) is the change in optical depth in direction n relative to the reference 
model. There is motivation for such a term to exist: if as is larger in some regions 
than others, the regions of higher as will form more early galaxies and we expect that 
they will reionize first. In the case where A is small, we may write: 

Q loc (n) « C, (1 + 2Xp • n) , (3) 

where the factor 

r lss d In a 8 

accounts for the increase of as and hence r in direction +p versus —p, and the factor 
r ri/ r iss acconts for the fact that the reionization surface is closer to us than the last 
scattering surface. Reionization occurs during the matter-dominated era when the 

3/2 

optical depth depends on the scale factor at reionization a ri via r oc a ri . Therefore 

Zrr r \ dlna r ; 
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From the WMAP 5-year parameter constraints we know that r = 0.087 ± 0.017 [2"7] 
and 3Tr r i/(2ri ss ) = 0.092 ± 0.021. One does not know d In a r ; /d In a$ , as this depends 
sensitively on astrophysics. If reionization is caused by massive haloes, so that the 
dependence of the reionization redshift on erg is dominated by changes in the halo mass 
function, then we have a r i oc a^ 1 because the halo mass function during the matter- 
dominated era is invariant under rescalings of a and as that hold asa constant. This 
case gives dlna r i/dlnr78 = — 1> i-e. higher as gives earlier reionization. This number 
will be increased somewhat (i.e. made closer to zero) if one takes into account that at 
earlier times more photons are required to reionizc the intergalactic medium because of 
recombinations. It will also be increased if higher as implies a larger clumping factor, 
or more minihalocs that act as a sink for ionizing photons. While the astrophysics 
is unclear at present, we are fortunate that the multiplying factor of 0.09 renders 
X relatively insensitive to the details: the d In a r j /d In eg = — 1 case would imply 
X w 0.91, whereas the extreme case of assuming reionization at an epoch independent 
of as would give X = 1. For practical purposes X = 1 with negligible correction 
(unless the anisotropy is measured at > 10c) and we will not consider it further. 

The power asymmetry analysis of Ref. [10] using the 3-year WMAP data and a 
range of scales 2 < £ < 40 found an asymmetry of A = 0.114 with the direction p given 
by (/, b) = (225°, —27°) and 99.1% significance. A more recent analysis by Hansen et 
al [IT] using the 5-year WMAP data and going out to ^ max = 600 finds preliminary 
evidence for a direction of (I, b) = (226° ± 10°, -17° ± 10°) at 99.6% significance. Very 
recently, Hoftuft et al [12] search the range of scales up to £ max — 64 and find a 3.3<r 
dipole, A = 0.072 ± 0.022 with p = (I, b) = (224°, -22°) ± 24° in the WMAP 5-year 
ILC map. In the 61 GHz band, they carried the analysis up to £ ma x = 80 and found 
A = 0.070 ± 0.019 (3.7cr) with p = (235°, -17°) ± 22°. 

An analysis by Donoghue & Donoghue [28] of the first acoustic peak finds a power 
spectrum asymmetry of r/ = 2A = 0.02 ± 0.02 at £ ~ 220 (assuming the direction p 
given by Hansen et al [9]), implying A < 0.03 at 2a. This is in tension with the low-£ 
results, and motivates further analysis to reduce the uncertainty on A. 

One way to reduce the error bars on the CMB method would be to extend it 
to higher £ where there are more modes available to constrain the hemispherical 
asymmetry. The number of CMB modes available is proportional to <^ ax and so 
we would expect the uncertainty in A to decrease as cx ^~ ax ; the Hoftuft et al 
uncertainties are in rough agreement with this expectation. Extrapolating this scaling, 
an analysis making use of the full WMAP data out to £ max ~ 500 should be able to 
reach a(A) ~ 0.003 if other sources of uncertainty (e.g. unmasked point sources) 
can be overcome. We believe such an analysis would be very helpful, especially to 
constrain the scale dependence of the asymmetry. However, in this paper we take an 
independent approach with similar sensitivity: we look for an asymmetry in large-scale 
structure. 

2.3. Effect on large scale structure 

The large-scale gradient in as also has an effect on the abundance of massive haloes and 
any objects that occupy them. For objects at a distance r, the observed 2-dimensional 
number density varies across the sky as: 




ri ss d In as 



r dh\N 



p ■ n. 



(6) 
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Our major problem here is determining d In N/d In cr 8 . This problem has been 
considered in the context of non-Gaussianity searches [351 EM OS] j which found 

SlniV , 

a In cr 8 

where <5 C = 1.69 and 6 is the bias. This formula assumed a universal mass function, 
as found in simulations, and assumed that the objects under study have a halo 
occupation distribution (HOD) that depends only on halo mass. The dependence 
on o~ 8 makes conceptual sense: very massive objects are highly biased b(M) ^> 1, and 
their abundance increases rapidly with erg; low-mass objects are antibiased b(M) < 1 
and are less abundant at high erg because they will have merged; and the overall 
abundance of dark matter particles with 6=1 does not depend on erg . If only recent 
major mergers are occupied, then there is an additional suppression in Eq. because 
massive haloes formed earlier if erg is increased. In the extended Press-Schechter (ePS) 
formalism, 6—1 should then be replaced by 6 — 1 — S^T 1 = b — 1.59 [31] ; this exponent 
has been confirmed by re-scaling of simulations [32l [33] which suggest that during the 
matter-dominated era the erg dependence for recent major mergers ranges between 
6 — 1.59 and 6 — 1.65 depending on the remnant mass and progenitor mass ratio |31j. 
In the case of the SDSS quasars it is not clear what is the correct way to populate 
haloes and so both limiting cases (6—1 and 6 — 1.6) must be considered. In practice 
the bias of the high-redshift quasars is so large (~ 10) that the difference between 
these cases is negligible. 

In practice the quasars occupy a range of redshifts and the observed dipole is a 
superposition: 5N/N = d ■ n, where 

d = 5 c p [ ^l[b(z)-l]P(z)dz, (8) 

with the replacement 6 — 1 — ► b — 1.6 in the recent major merger case. Here P(z) is 
the redshift distribution of the quasars, normalized to unity: J P{z) dz = 1. 



3. Application to SDSS quasars 

We now turn to empirical determination of Vug using the above formalism. 



3.1. Choice of dataset 

The high-redshift spectroscopic quasars from SDSS were chosen for this study because: 

(i) They are very distant, which means that the same spatial gradient in erg translates 
into a large absolute difference Aug across the sky. 

(ii) They have wide angular coverage, which mimimizes large-scale structure noise 
and provides the leverage in all three Cartesian coordinates necessary to measure 
a dipole. 

(iii) They are highly biased, so that a small change in erg is amplified into a much 
larger change in the number density of quasars. 

(iv) They have very good rejection of Galactic objects (e.g. halo stars), which 
could otherwise produce spurious large angular scale modulation and potentially 
correlate with foreground residuals in the CMB maps. 

(v) The relative photometric calibration of SDSS is understood at the percent level 
(although this degrades to a few percent in the u band). 
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(vi) The number density of ~ 1 deg~ 2 , while not impressive, is acceptable given 
the above advantages. Given the choice between high number density of nearby 
objects with a weak dependence on er 8 , or a low number density of distant objects 
that depends strongly on er s , our preference is for the latter since Nature does 
much of the work to suppress systematics. 

There are other data sets that could have been used instead to trace large-scale 
structure at moderate to high redshifts. Possibilities include: 

(i) Luminous red galaxies: These are abundant but current large-scale maps go out 
to only z < 0.7 [34j . or 18% of the distance to the last-scattering surface, and the 
typical bias is small, b ~ 1.8. 

(ii) Photometrically selected quasars: These have a higher number density (~50 
deg -2 ), but at present the UV-based selection algorithms restrict this sample 
to z < 2.5 |35| . The bias is also much lower at these redshifts, b ~2.3-2.8 [55] , 

(iii) Hard X-ray background: The 2-8 keV HEAO map has excellent sky coverage and 
perhaps the best-constrained dipole [37], but the redshift distribution is poorly 
known and the estimated bias is small and may be consistent with 6 = 1 [38] , so 
it is not known how sensitive this is to Vers. 

(iv) Radio sources: The NRAO VLA Sky Survey [33] sample has been used to probe 
the moderate redshift range z ~ 1, but its redshift distribution is still under debate 
[401 1411 136j . More importantly the maps exhibit declination-dependent striping 
which is almost certainly a systematic artifact and precludes determination of the 
dipole P] . 

(v) Lyman-a forest: A different approach would be to measure a fluctuation 
amplitude from the Lyman-a forest at the same redshift but with sightlines in 
different parts of the sky. The differential nature of the measurement would 
suppress many of the astrophysical systematics in e.g. P(k) measurement from 
the Lyman-a forest. This approach was not taken here because the analysis would 
be much more complicated, but given the large number of observed modes in the 
Lyman-a forest it should be considered in the future. 

3.2. Data description 

The high-redshift quasars used here are obtained from the SDSS. The SDSS drift- 
scans the sky in five bands (ugriz) [43] under photometric conditions [44] 05] using a 
2.5-meter optical telescope [46 with 3 degree field of view camera [47] located in New 
Mexico, USA [Uj. The photometric and astrometric calibration of the SDSS and the 
quality assessment pipeline are described by Refs. [1H1 HS1 H33 US [52] • Bright galaxies 
[53] , luminous red galaxies (LRGs) [54] , and quasar candidates [55] are selected from 
the SDSS imaging data for spectroscopic follow-up [56] , 

The spectroscopic quasar sample used here was based on colour selection criteria 
described in Richards et al [55]. The algorithm selects quasar candidates down to 
a limiting magnitude of i = 19.1 or 20.2 (for high-z quasars); a discussion of its 
completeness and efficiency can be found in Richards et al [55]. A statistical sample 
of confirmed quasars was generated by Shen et al [57] , who removed from the sample 
regions with early versions of the target selection algorithm. Shen et al imposed a 
further cut on the quality of the imaging data used for target selection, dividing it into 
"good" and "bad" fields. As a test for systematics, their quasar clustering analysis 
was repeated with both the "good" fields only, and with "all" fields. 
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The Shen et al [57| sample has a minimum redshift z > 2.9, due to the need in 
target selection to avoid Galactic late A/early F stars which have similar broadband 
colours to quasars at z ~ 2.7, and most objects are at z < 4.5 due to the declining 
luminosity function and increasing luminosity distance at high redshift. There are 
4426 quasars in the "all" sample (4041 deg 2 ) and 3846 in "good" (3506 deg 2 ). The 
quasars undersample the density field at all scales, so that their power spectrum is 
dominated by Poisson noise. 

Shen et al specify a mask, which we have imported into Mangle [53 \59\ [60] . 
Mangle is a suite of computer programmes widely used in the large-scale structure 
community to manipulate complex masks, which it typically represents as spherical 
polygons (i.e. regions bounded by circular arcs). Mangle can carry out simple tasks 
on these masks such as unions and spherical harmonic transform computation; most 
importantly for us, it can generate catalogues of random points from within a given 
mask. 

For this analysis, we have further split the quasars into a 2.9 < z < 3.5 sample, 
and a 3.5 < z < 4.5 sample. 



3.3. Dipole anisotropy estimator 

In general the quasar density in a given direction can be written as: 

^(») = d-n + ^2kiti(h) + C. (9) 

i 

Here N is the mean number of quasars per steradian, SN is its fluctuation, d is the 
dipole (which we wish to measure), tj(fi) are possible systematics templates in the 
quasar maps (e.g. the extinction map), ki are the sensitivities of the quasar densities 
to these systematics, and C is a mean offset over which we must marginalize since 
the true mean number density of quasars is not known. This equation can be simply 
written as: 

^(n)=*.T(n), (10) 

where we have defined the vectors x — (d, ki, C) and T(h) = (n,ti(n), 1). If Poisson 
noise dominates the uncertainty in SN/N, so that the inverse variance per steradian 
is N, then the best linear unbiased estimator for x is obtained via: 

x = F- 1 g, (11) 

where 

r Ti(n)SN(h)d 2 h (12) 
and the Fisher matrix is 

F tj = N I r j; (n)r,-(n)d 2 n. (13) 



The covariance matrix of this estimator is Cov(x) = F _1 . 

The above equations are integrals over the survey area, but quasars are discrete 
objects and so it is easiest to replace the above results with summations over the 
quasars and random catalogues. We can do this by writing: 

D N *R 
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where rio are the positions of the quasars in the real data, are the positions of 
randomly generated points in the survey mask, and iVrj and are the number of 
quasars (data) and random points respectively. The Fisher matrix is: 



We have generated random catalogues containing TVr = 10 6 points using Mangle's 
ransack task [60], which is more than sufficient as their residual noise is negligible 
(< 1%) compared to the actual data. 

3-4- Measured dipole anisotropy 

We have applied the estimator Eq. (fTT|) to the high- 2 SDSS quasar data. Four versions 
of the fit were done, which differed by using the 'all" or "good" field mask, and by either 
including or not including the dust template. The dust template used is the predicted 
reddening map of E(B — V) by Schlegel et al 61j as a systematics template, which is 
appropriate if the extinction correction is not perfect. Errors correlated with the map 
could arise due to errors in the assumed extinction law A\/E(B — V), miscalibration of 
the E(B — V) scale, or unusual quasar colours affecting the extinction correction (e.g. 
if a bright emission line lies at the red end of the bandpass then the band-averaged 
correction may overcorrect for the extinction). The template amplitude kE(B-v) could 
plausibly be positive or negative depending on the amount of over- or undercorrection 
in different bands. Because most of the data lie in the Northern Galactic cap, there is a 
strong degeneracy between the component of the dipole perpendicular to the Galactic 
plane and the extinction correction, Cott(kE(B-V)t d z ) = 0.65: either increasing d z 
or decreasing k^m-V) would have the effect of enhancing the number density of 
objects near the NGP. This degeneracy could be broken in the future with more 
data in the southern patch. For the "all" mask, we find template amplitudes of 
k,E(B-v) = —1-6 ± 1.2 (low-z) and k.E{B-v) = —3.5 ± 1.6 (high-z); these become 
— 1.0 ± 1.3 and —3.6 ± 1.7 for the "good" mask. We regard the fits with the dust 
template as our primary result because correlation with dust is marginally detected 
at > 2d in one of the redshift slices. 

The resulting dipole moments in the Galactic coordinate system are shown in 
Table [TJ To be explicit, we have placed the x direction toward the Galactic centre, 
the y direction in the plane at I = 90°, and the z direction toward the North Galactic 
Pole (NGP). Using the "all" mask and no dust template, one finds an intriguing 2.9<r 
hint of d z > in the low- z quasars. However the fact that the preferred direction is 
very near the NGP (b = 80°) is suspicious, and indeed this signal goes away if one 
restricts to "good" fields (significance drops to 1.8ct), uses the dust template (1.4er), 
or both (0.8ct). 

In all cases the uncertainty in the dipole determination is at least several percent, 
so no correction for our own peculiar velocity v/c = 0.0012 [2 is necessary. 

In addition to the Galactic dust, a second concern would be large-angle calibration 
errors, which can involve either "grey" errors (in which all five bands move together) 
or colour calibration errors. We may do a simple assessment of the grey calibration 
error by comparing the target selection photometry to the ubercalibrated photomety 
[50] obtained from the Shen et al [57, catalogue, in i band as used for the flux cut. 
The ubercalibrated photometry is estimated to be accurate to 1%. We compare the 




(15) 
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Table 1. The observed dipole anisotropies, in Galactic Cartesian coordinates, 
for the various quasar samples and for the dust template either included (Yes) or 



not (No). 


Sample 


Dust? 


I0 2 d x 


lO 2 ^ 


10 2 d z 




2.9-3.5, all 


No 


- 1.3 ±5.1 


2.6 ±6.5 


16.3 ± 


5.6 


3.5-4.5, all 


No 


- 0.5 ±6.7 


- 4.2 ±8.4 


- 1.9 ± 


7.3 


2.9-3.5, all 


Yes 


1.1 ± 5.4 


- 0.7 ±6.9 


10.0 ± 


7.1 


3.5-4.5, all 


Yes 


4.8 ± 7.1 


-11.7 ±9.0 


-15.9 ± 


9.6 


2.9-3.5, good 


No 


- 2.3 ±5.5 


2.8 ± 6.9 


10.8 ± 


6.0 


3.5-4.5, good 


No 


3.2 ± 7.1 


- 2.7 ±9.0 


- 7.4 ± 


6.8 


2.9-3.5, good 


Yes 


3.9 ±5.8 


0.6 ± 7.5 


6.6 ± 


7.9 


3.5—4.5, good 


Yes 


8.7 ±7.6 


-10.5 ± 9.6 


-21.9 ± 


10.2 



magnitude difference Ai = itargct — iubcrcai for each quasar, and do an unweighted 
least-squares fit 

Ai = Aio + m ■ n, (16) 

where m is the calibration dipole and n is the position of that quasar. We remove 
two (possibly variable) objects from the fit that were undetected in the ubercalibrated 
photometry. The resulting dipole is 10 3 m = (0 ± 3, — 10 ± 3, — 3 ± 3). We have also 
tried clipping the distribution of Ai at the 1st and 99th percentiles (i.e. quasars below 
the 1st percentile were moved to the 1st percentile value before fitting). The resulting 
dipole is 10 3 m = (0.6 ± 0.8, -3.5 ± 0.9, 0.5 ± 0.8). The implied dipole due to changes 
in the i band calibration would be 1.8 times this, since the slope of the cumulative 
quasar luminosity function for this redshift range and the i < 20.2 flux cut is ~ 1.8 
(see e.g. Fig. 13 of Ref. [62]). Even the 1% dipole obtained from the first fit would 
lead to a calibration-induced dipole of 1.8% in the density of quasars, which is much 
less than our 5-10% statistical errors. (But do note that both calibration dipoles are 
inconsistent with zero, giving x 2 for 3 degrees of freedom of 12 and 19, respectively.) 

Colour calibration is potentially nastier since quasars are selected by cuts that run 
through the quasar locus, and hence could result in completeness variations if there is 
a relative colour calibration error. Our best test of the colour-dependent calibration 
will be the agreement (or lack thereof) of the dipoles inferred from different redshift 
ranges where the colour selection is very different. However, we can also test for this 
by computing the dipole of the reddening-corrected colours, i.e. doing an unweighted 
least-squares fit 

u-g = (u-g) + m-n, (17) 

where (u — g)o and the vector m represent 4 free parameters. We use all of the quasars 
in the 2.9 < z < 4.5 range for this test. Since some quasars are undetected in u, we 
clip the distribution at the 5th and 95th percentiles in u — g before performing the fit 
(i.e. quasars below the 5th percentile were moved to the 5th percentile value before 
fitting). One can then test for whether m ^ 0, which would possibly indicate a dipole 
variation in colour calibration across the sky. We report the % 2 value, m ■ Cov^^m, 
and the probability to exceed P(> \ 2 )- Similar tests can be done for the other colours, 
and for the extinction-corrected apparent magnitudes of the quasars. The magnitudes 
and colours used are those used for the target selection (TARGET PSF photometry), 
which were obtained via matching to the SDSS DR5 quasar catalogue [B3]. The results 
are shown in Table [2j As one can see, the x 2 values for these fits are all acceptable, 
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Table 2. The \ 2 values and probabilities to exceed for searches for dipolc 
gradients in the quasar colour. 



Quantity 


X 2 (m) 


P(>X 2 ) 


u - g 


3.05 


0.38 


g-r 


4.61 


0.20 


r — i 


9.10 


0.028 


i — z 


2.84 


0.42 


u — r 


2.22 


0.53 


9-i 


4.53 


0.21 


r — z 


6.29 


0.10 


u — i 


1.58 


0.66 


g-z 


4.92 


0.18 


u — z 


1.36 


0.72 



with one possible exception (r — i, p = 0.028). Given that we searched 10 directions 
in colour space, the presence of this one anomalously high \ 2 value is not a concern. 

3. 5. Quasar bias and final results 

In order to convert the measured dipole anisotropy into a constraint on the 
perturbation amplitude asymmetry A using Eq. ([5]), we need to know the bias of 
the quasars. If one assumes a cosmological model, one may obtain this by fitting 
to the correlation function. Here we use the WMAP5+BAO cosmology, based 
on combining the WMAP 5-year power spectrum [64] with the SDSS+2dF baryon 
acoustic oscillation (BAO) constraint [65]; this combination yields as = 0.807 ±0.027 

If one has a measured correlation length ro and an assumed cosmology, the 
bias can be obtained by the formula b = [S.m^o)]^ 1 ^ 2 , where £ m (r) is the matter 
correlation function. We have used the linear correlation function using the Eisenstein 
& Hu transfer function |66j since in the regime of interest the nonlinear corrections as 
estimated by the Smith et al 67J prescription are < 1%, i.e. negligible. The correlation 
lengths that Shen et al [57] measured from the "good" fields are 17.91 ± 1.50 hr 1 Mpc 
(2.9 < z < 3.5) and 25.22 ± 2.50 h^ 1 Mpc (z > 3.5). These correlation lengths imply 
b = 8.9 ± 0.7 (z = 3.2) and b = 15.3 ±1.7 (z = 4.0), where the uncertainty is 
simply propagated from ro. There is some additional uncertainty in these numbers 
from nonlinear biasing, from the cosmological parameters, and because the correlation 
length does not appear to be constant with redshift. We note that the one-halo term 
in the correlation function cannot contaminate these measurements because Shen et 
al only fit projected separations r p > 4/i~ 1 Mpc. Integrating over the redshift bin 
assuming a constant correlation length r , we turn Eq. ([5]) into: 

(5.8 ± 0.6)p low— z slice , . 

\ (11.5±1.4)p high- z slice ' [ ' 

Using the central value of the conversion from the observed quasar dipole into p, 
and the conservative (good + dust marginalization) dipole estimates, we obtain 10 2 p = 
(0.7 ± 1.0, 0.1 ± 1.3, 1.1 ± 1.4) for the low-z slice and (0.8 ± 0.7, -0.9 ± 0.8, -1.9 ± 0.9) 
for the high-z slice. Combining the two results with the usual weighting by their 
inverse-covariance matrix gives 

10 2 p = (0.73, -0.61, -1.01), (19) 
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with covariance: 

/ 0.306 -0.225 -0.122 

10 4 Cov (p) = -0.225 0.497 0.227 ] . (20) 

\ -0.122 0.227 0.560 

In the optimal determination of the central value fEq. [T9|) the low-redshift slice receives 
30% of the weight and the high-redshift slice receives 70%, i.e. 

P = Wl _ z pio- z + Whi-zPhi-z (21) 

with w\o- z — 0.30 and Whi- Z = 0.70. 

The difference between the inferred dipoles from the two samples is pi - z — Pm-z = 
(—0.1, 1.0, 3.1) x 10~ 2 ; based on the covariance matrix for the difference, the x 2 1S 3.03 
(3 degrees of freedom; P = 0.39) so we conclude the dipoles measured from the two 
redshift slices are consistent. There is no detection of any global dipole (x 2 = 2.80; 
P = 0.42). 

If one fixes p to be in the specific direction (l,b) = (225°,— 27°) identified by 
Eriksen et al 10J, then the magnitude A of p is found to be 10 2 A = -0.18 ± 0.44, 
i.e. -0.0105 < A < 0.0070 at 95% confidence (-0.0132 < A < 0.0097 at 99%). 
A directon-independent upper limit can be set by a Bayesian analysis analogous to 
Eriksen et al : by placing a uniform prior on the magnitude A and direction p, we 
can construct a marginalized likelihood function for A: 



C{A) oc J exp 



" 2 (P ~ Pbcst ) • COV 1 (p - Pbcst ) 



d 2 p, (22) 



where p = Ap. We find that 95% of the Bayesian posterior distribution is at A < 0.019 
and 99% at A < 0.026[jJ 

In reality the above computations should take into accounte the uncertainty in 
the estimator response, Eq. (|18[) . which derives from the uncertainty in the quasar 
bias. In particular it is possible to circumvent the above limits on p if the quasar bias 
is smaller than central values. For our chosen weights, Eq. (f2Tj) . the response of our 
estimator is: 

^estimated = (1-00 ± 0.09)p actua l, (23) 

i.e. there is a 9% (lc) calibration uncertainty. We can include this uncertainty in a 
Bayesian analysis by incorporating a calibration factor x with a prior at 1.00 ± 0.09. 
For the fixed-direction case (i.e. where p is required to point in the Eriksen et al 
direction), one may define a margimalizcd likelihood function for A: 

-(x^+0.0018) 2 /2(0.0044) 2 -(x-l) 2 /2(0.09) 



Anarg(A) CX / e -W+™^*r /^.UU44 r ^-1)-/^-^)- dx . ( 2 4) 

This gives -0.0108 < A < 0.0071 (95%) or -0.0138 < A < 0.0100 (99%). For the 
direction-independent case, 

Anaxg(^) = / A m - marg (xv4)e-(^ 1 ) 2 / 2 (°- 09 ) 2 d x , (25) 



where the un-marginalized function £un-marg (A) is given by Eq. (f2"2")) . Here we 
find that 95% of the Bayesian posterior distribution is at A < 0.020 and 99% is 
at A < 0.027. These are negligible changes from the un-marginalized case. 

|| Since A is the magnitude of p, in our direction-independent constraint it is required to be non- 
negative. In the fixed-direction constraint where p is fixed to the Eriksen et al direction, we allow 
for the possibility of negative A, i.e. a dipole opposite to that seen in the CMB. 
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We have also repeated the above analysis for the "all" fields. In this case, the 
correlation length measured by Shen et al 57J is shorter, so the biases are less: 
b = 8.0 ± 0.9 at z = 3.2 and b — 13.5 ± 1.6 at z = 4.0. The dipole response is 
now 

, J (5.1±0.7)p low— z slice , , 

\ (10.0 ± 1.3)p high-z slice ' ' ' 

The best-fit dipole and its covariance are: 

10 2 p = (0.40, -0.87, -0.52) (27) 

and 

/ 0.351 -0.259 -0.141 \ 
10 4 Cov (p) = -0.259 0.571 0.261 . (28) 
\ -0.141 0.261 0.643 / 

The weights in this case are wi D _ z = 0.30 and u>hi- z = 0.70. Again the low-z and 
high-z slices are consistent (\ 2 = 3.11, P = 0.37). There is also no detection of p: the 
X 2 for zero dipole is 1.34 (P = 0.72). The calibration uncertainty is now 10%. After 
marginalizing over this, the fixed-direction constraint on A is —0.0073 < A < 0.0120 
(95%) or —0.0105 < A < 0.0153 (99%), and the direction-independent constraint is 
A < 0.017 (95%) or A < 0.025 (99%). 

3. 6. Scale of constraint 

Up until now we have assumed that the dipole anisotropy p is scale- invariant. However 
inflationary models usually predict slight deviations from scale invariance because as 
one approaches the end of inflation each e-fold is slightly different from the previous 
one. In order to constrain such models, we need to know the effective scale at which the 
quasars constrain p. Since the fluctuation amplitude is typically a smoothly varying 
function of the number of e-folds or equivalently of Ink, only a rough estimate of fc e ff 
is required. 

We now obtain an estimate of the scale probed by the quasars, i.e. at what scale 
fc e ff we have actually constrained the large-scale gradient of the fluctuation amplitude. 
The abundance of very massive haloes (i.e. well above the nonlinear mass) is roughly 
determined by the top-hat variance cr(R), where M = 47rp m oi? 3 /3 is the halo mass, 
p m Q is the present-day matter density, and R is the comoving radius of the top-hat 
filter. This is: 

a 2 {R) = J A 2 (k)W 2 {kR)d\nk, (29) 

where the window function is W(x) — 3j\(x)/x. Differentiating gives 

/A 2 (fc)^ 2 (fci?)JlnA(fc)dlnfc 
61na(R)= J A2{k)W 2 (kR) dlnk • ( 30 ) 

This can be approximated as 5 In o~{R) = iJln A.(k e g), where: 

_ J A 2 {k)W 2 {kR)\nkd\nk 
eff - J A2 (k)W 2 (kR) dlnfc ' (31) 

For the quasars, the minimum mass (which for a steep mass function can be taken as a 
typical mass) is M ~ (2-6) X 10 12 /i _1 M Q [57], which gives k eS = (1.3- 1.8)/iMpc _1 . 
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This range spans only 0.3 e-folds and in what follows we will thus take the central 
value, k e s — 1.5ft Mpc -1 . 

One can make a similar estimate for the Eriksen et al |10j ; we find in | Appendix A| 
that fc e ff = 0.0033/i Mpc -1 . The Eriksen et al constraint and our constraint are thus 
separated by 6.1 e-folds. This strongly constrains inflationary explanations for the 
power asymmetry: the asymmetry amplitude A must decrease by a factor of at least 
a few (~ 7% to <~ 1%) during only 6.1 e-folds of expansion. This contrasts with the 
usual expectation from inflationary models of nearly scale-invariant spectra. One could 
imagine features in the inflationary potential that strongly break scale invariance; 
however there is no motivation to place such features in precisely this range of In fc, 
and moreover such a model would have to avoid "breaking" the successes of standard 
inflation in predicting the CMB acoustic peak regime and large-scale galaxy clustering. 
A detailed study of such constraints is presented in a companion paper |69j . 

4. Discussion 

One of the intriguing products of the WMAP mission was the detection of large- 
scale anomalies in the CMB. If these anomalies are "real" (as opposed to statistical 
fluctuations or systematics) then their implications are revolutionary. The anomalies 
are thus in need of confirmation, both by future CMB data that could detect them 
at higher signal-to-noise ratio (e.g. Planck) and just as importantly by independent 
probes of the same underlying physics. In particular, the dipolar asymmetry in the 
CMB power spectrum reported by Eriksen et al [TU], if interpreted as a spatial 
variation of the amplitude of primordial fluctuations, makes predictions not just for 
the CMB but also for large-scale structure: if there is truly a large-scale gradient 
Vers across the Universe, then there shuld be a corresponding gradient in the number 
density of highly biased objects. We have searched for such a gradient using the SDSS 
z > 2.9 quasars and found a null result. We find that the dipole p = ri ss Vlnc78 
is no more than 0.027 (99% posterior probability), i.e. any smooth gradient in the 
amplitude of fluctuations is no more than 2.7 percent per present-day horizon radius. 
Tighter results are obtained if one forces the gradient to be in the same direction as 
that reported by Eriksen et al |10j . 

This result rules out the simple curvaton-gradient model |20j in which the power 
asymmetry in the CMB comes from a spatial gradient in the curvaton field of a two- 
field inflation model. Such a model requires p w 0.11 to reproduce the reported 
±11% variation in the CMB amplitude across the sky rjUJ. This does not necessarily 
contradict the Eriksen et al [10] result: there may be more complicated models with 
a scale-dependent gradient in small-scale power, i.e. where the spatial gradient of the 
local lnP(fc) depends on k. Such models, including some that predict much smaller 
dipoles in the quasar distribution and the degree-scale CMB fluctuation amplitudes, 
are explored in the companion paper by Erickcek et al [69 . 
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Appendix A. Effective scale of CMB power asymmetry 

Eriksen et al [lOj measured a dipole asymmetry in the CMB power spectrum in 
the range of multipoles 2 < £ < £ max ~ 40. The upper limit i max was set by the 
requirement to be oversampled on the Healpix resolution 4 (A s id c = 16) pixelization 
scheme [55]. Our goal here is to estimate the effective wavenumber of the Eriksen et 
al measurement. 

In the Sachs-Wolfe regime, and neglecting the small contribution from the 
integrated Sachs- Wolfe effect, the temperature perturbation is given by Q(n) — 
— C( r iss^)/5 where £ is the primordial curvature perturbation. Therefore the CMB 
power spectrum at multipole I is obtained by 
Ait f 

Ce =25j 3i(kn ss )Al(k)dlnk, (A.l) 

where A^(fc) is the primordial curvature power spectrum. [See, e.g. Eqs. (5.27) and 
(5.39) of Ref. [70] ■] For smooth power spectra and f > 1, we can use the Wentzel- 
Kramers-Brillouin approximation for the spherical Bessel function, 

j £ (x) -> x^/^x 2 -i 2 )- 1 / 4 cos ip(x) (A.2) 

for x > £, where <f{x) is a rapidly oscillating phase. (For x < I the spherical Bessel 
function goes rapidly to zero.) With this approximation, and taking (cos 2 f{x)) — -> h, 
Cg simplifies to 

G ^Mp\ V ~ 1{y2 - ^^(kWnk, (A.3) 
where y — xjl — kr\ ss /£. The fractional dipole asymmetry SlnCe of Ci is given by: 

/y-V-ir 1/2 A 2 (£0<51nA 2 (fc)dlnfc 
6lU Ce * jVi(y 2 -l)-V2 A 2 (fc)dlnfc ■ ( A ' 4) 
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This equation can be simplified for the case of a near scale-invariant reference 



spectrum, A?(fc) ^constant: 



"" c "* c»-v-D-^i°» ' <A ' 5) 



The integral in the denominator evaluates to 1 (with the substitution y = sec/3) 



In order to make further progress, we assume that 5 In A? (A;) varies smoothly with 



Ink, as typical during inflation models, i.e. we assume: 

S\nA 2 c (k) = B x + S 2 ln(fcr lss ) = Si + S 2 hil + S 2 lny. (A.6) 
Only the In y part of the Ci integral is nontrivial, 

/>oo 

6\nC e ^B 1 +B 2 ln£ + B 2 / jtV - l)" 1 / 2 lny dlny 

Jy=l 

= S x + S 2 ln^ + S 2 (l -ln2). (A.7) 

{The integral over y can be solved by the substitution y = (1 — p 2 )" 1 / 2 , which turns 
it into — i J* [ln(l — p) + ln(l + p)]dp. The In functions are then integrated by parts.} 
The Eriksen et al [10] dipole asymmetry measurement fits a single amplitude for 
the power asymmetry across their entire range of £. In the limit where all multipoles up 
to f m ax are signal-dominated and those above £ max are noise-dominated, the Eriksen 
et al asymmetry parameter A will receive equal weight from all modes, i.e. 

A = -(SlnC e ) modes » i [Si + S 2 (ln£) modcs + S 2 (l - In 2)] . (A.8) 

(The factor of 1/2 comes from the fact that Eriksen et al measured an amplitude 
asymmetry rather than a power asymmetry.) Since the number of modes in any 
region of the sky and any £ range is proportional to £d£, we have, for ^ max 3> 1, 

(ln^) mod e 8 =^ =^£ max --. (A.9) 

J—£d£ 2 

(The integral in the numerator is evaluated using repeated integration by parts.) Thus: 

(A.IO) 



A = -(5 In Ci) modes ~ ^ 



Lh K/->> ! In£ max + i-ln2 



This is the same as the amplitude asymmetry ^5 In A 2 (fc c ff) at the effective scale: 

fc off = e -J^, (A.11) 
2ri ss 

where c = 2.718... is the base of the natural logarithm. 

For the specific case of Eriksen et al , n ss = 1.0 x 10 4 /i _1 Mpc and ^ max = 40, 
so fceff = 0.0033/j Mpc" 1 . Similarly, for Hoftuft et al [12], k eS = 0.0053/i Mpc" 1 
(4iax = 64) or fees = 0.0066/i Mpc" 1 (^ max = 80). 
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